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ABSTRACT 

We investigate the ratio between the half-mass radii of Galactic globular clusters 
and their Jacobi radii rj given by the potential of the Milky Way and show that 
clusters with galactocentric distances Rgc > 8 kpc fall into two distinct groups: one 
group of compact, tidally-underfilling clusters with ri^/rj < 0.05 and another group of 
tidally filling clusters which have 0.1 < ru/rj < 0.3. We find no correlation between 
the membership of a particular cluster to one of these groups and its membership 
in the old or younger halo population. Based on the relaxation times and orbits of 
the clusters, we argue that compact clusters and most clusters in the inner Milky 
Way were born compact with half-mass radii r/j < 1 pc. Some of the tidally-filling 
clusters might have formed compact as well, but the majority likely formed with large 
half-mass radii. Galactic globular clusters therefore show a similar dichotomy as was 
recently found for globular clusters in dwarf galaxies and for young star clusters in the 
Milky Way. It seems likely that some of the tidally-filling clusters are evolving along 
the main sequence line of clusters recently discovered by Kiipper et al. (2008) and are 
in the process of dissolution. 
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1 INTRODUCTION 

It is well known that most, if not all, stars form in star clus- 
ters. Star clusters are therefore important probes o f the star 
formation process l|Kroupall200"5l : |Parmentieij|2009l ). This is 
especially the case for globular clusters which, due to their 
large ages and low metallicities, are relics of star formation 
processes in the early universe. 

Observations show that in nearby galaxies star clus- 
ters form compact, with half-mass radii of < 1 pc, at 
the c entres of giant molecular cloud cores (|Lada fc Ladal 
They then undergo a significant expansion as a re- 
sult of gas expulsion driven by stellar winds and U V ra- 
diatio n from bright stars an d sup ernovac explosions fHills! 
198C ; [Bastian fc Goodwiry l2006l : jBaumgardt & Kro upa 
20071) and later also by mass loss from stellar evolu- 



tion llChernoff fc Weinbe"r3ll990l : iFukushige fc HeggidflQQSI : 
IVesperini fc Zepj|2003l )." As a result, the radii of star clusters 
show a steady incr ease with cluster age within the first 20 
Myrs (jBastian et al . 2008; Pfalznor 2009). 

Globular clusters are dynamically evolved systems, 
so in addition to the early evolution due to gas ex- 
pulsion and stellar evolution, cluster radii are also sub- 
ject to dynamical cluster evolution due to two-body re- 



laxation. As a result of mass segregation and core col- 
lapse, the core radius shrinks while the half mass ra- 
dius stays roughly constant before core collapse. After 
core collapse, stellar binaries provide a ce ntral heat source 
and the cluste r expands self-similarly (iGoodm anI Il984|; 
McMillan et~al] Il990l: iGao et al.l Il99ll: iGiers z & Hcgdi 
19941 : iBaumgardt. Hut fc Heggid I2OO2I : iHeggie et al.i 120061 )? 
This process continues until the cluster runs into the ex- 
ternal tidal field, at which point cluster expansion is bal- 
anced by the loss of outer stars over the tidal boundary. 
As a result, the ratio of half-mass radius to Jacobi ra- 
diu^ evolves along a common sequence which, at least for 
single-mass clusters, only depends on current cluster mass 
jKupper. Kroupa fc Baumgardtll2008h . The ratio of the half- 
mass to Jacobi radius for individual clusters therefore con- 
tains important information on the dynamical state of a star 
cluster. 

In the current paper we examine the distribution of half- 



^ Throughout the paper we will denote the distance from the 
centre of a star cluster to the first Lagrangian point as the Jacobi 
radius rj, while the term tidal radius rt refers to the limiting 
radius of King (1962) or King (1966) models. 
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mass to Jacobi radii of Galactic globular clusters in order to 
better understand the formation and evolution of Galactic 
globular clusters. The paper is organised as follows: In Sec. [2] 
we present the data and show that outer globular clusters 
fall into two distinct groups. In Sec. [3] we discuss the relation 
of these groups to different subsamples of Galactic globular 
clusters like old and younger halo clusters or globular clus- 
ters believed to be accreted from dwarf galaxies and use 
additional information like half-mass relaxation times and 
cluster orbits to determine the degree of dynamical cluster 
evolution. In Sec. |4]we finally draw our conclusions. 



2 GLOBULAR CLUSTER DATA 

In order to investigate the degree of tidal filling, one has to 
obtain an estimate of the Jacobi radius of a cluster. Most in- 
vestigations so far used the tidal radius rt obtained by fitting 
the observed sur face density profile with an empirical profile 
like lKinid (|l962l ) or a theoretical one like lKind l|l966h as an 
estimate of the Jacobi radius. Surface density data for most 
globular clusters is however either not available near the Ja- 
cobi radius or becomes unreliable in the outer parts due to 
the low number of cluster stars and the uncertain density 
of background stars. In addition, due to the high number 
of stars, surface densities can be much better determined 
in the inner cluster parts, so that the published tidal radii 
rt of globular clusters are determined more by the density 
profile inside a few half-mass radii and mi ght not reflect the 
true tidal radius of a cluster. Baumgard t et al.l (|2009l ) for 
example found that in case of NGC 2419, the best-fitting 
King model has a nominal tidal radius of 150 pc, while they 
estimated that the Jacobi radius r j of the cluster is around 
800 pc. 

An additional problem of fitting King models to deter- 
mine Jacobi radii is depicted in Fig. [T] This figure shows 
the ratio between the proje c ted h alf-m ass ra d ius a nd the 
tidal radius rhp/n for iKind l| 19621 ) and iKinel (|l966l ) mod- 
els of various concentrations c — log-^^grt/rc■ It can be seen 
that with both families of mo dels only a limited range of 
Thp/rt values can be reached. iKingj l|l962l ) models are re- 
stricted to values between 0.028 ^ r^p/rt ^ 0.29 while 
iKind l| 19661 ) models are restricted to the even smaller range 
0.074 ^ Thp/rt ^ 0.23. A typical globular cluster with a 
mass of Mc = 2 • 10^ Mq and a projected half-mass radius of 
rtip = 3 pc has rhp/rj = 0.029 at a galactocentric distance 
of Rgc = 10 kpc and rup/rj — 0.018 at Rac = 20 kpc. 
Hence King models cannot accurately describe the density 
profile of such a cluster and one would have to use models 
allowin g for m ore extended envelopes like those described by 
IWilsonI (I1975I ) to fit the outer surface density profile . This 
was also noted bv lMcLaughlin fc van der Marell (|2005l ). who 
found that Wilson models provide equally good or signifi- 
cantly better fits than King (1966) models for about 90% 
of their sample of young massive clusters and old globular 
clusters. 

In the present paper we use the ratio of the half-mass 
radius to the Jacobi radius r j as a measure for the degree of 
tidal filling that a globular cluster experiences in the tidal 
field of the Galaxy. The Jacobi radius can be determined if 
the cluster mass, galactocentric distance and the underlying 
Galactic potential is known. For typical globular clusters. 




M^"=2- IQS M^; "R,,;,="20 kpc 
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Figure 1. Ratio of projected half-mass radius r^p to tidal radius 
rt as a function of the concentration index c = logj^Q rt/vc for King 
(1962) and King (1966) models. The dotted lines mark expected 
ratios of r^^p/rj for a typical globular cluster with Mc = 2 ■ 
10^ Mq and r^p = 3 pc at Rqc = 10 kpc and Rac = 20 
kpc galactocentric distance. The r^^^lrj ratio of such a cluster 
cannot be reproduced by any King model, showing that fitting 
King models to observed density profiles in order to derive the 
Jacobi radius can lead to a significant bias towards a too small 
value. 

masses and galactocentric distances should have typical er- 
rors of less than 20%, hence Jacobi radii can for many clus- 
ters be determined with higher accuracy than their tidal 
radii rt. The main disadvantage is that the Jacobi radius 
varies along the orbit of a cluster, so for highly eccentric or- 
bits, the current radius might not be a good measure for the 
average Jacobi radius which a cluster experiences and which 
determines its mass loss rate. We will discuss the influence 
of eccentric orbits on our results further below. 

In order to calculate the ratio r^/rj of Galactic globu- 
lar clusters, we calculated 3D half-mass radii rh from the 
projected half-light radii r^p under the assumption that 
mass follows light and by assuming rn = 1.33r^p. T his rela- 
tion is correct to within 5% for most iKind l|l962l ) or iKind 
(|l966t ) density proflles. We note that half-light radii can 
be different from half-mass radii for highly evolved clusters 
that have undergone core -collapse (|Baumgardt fc Makind 
l2003l : [Balbinot et al.ll2009l ). However, correcting for this ef- 
fect would require detailed observational data for each clus- 
ter which is not availa ble at t he mo ment. Jacobi radii were 
calculated according to lKind (Il962 l) (see also Innanen, Har- 
ris & Webbink 1983 who added a factor of 2/3 to correct for 
the elongation in the direction along the line connecting the 
Lagrangian points): 

Here is the mass of the cluster, Vg the circular velocity 
of the galaxy and Rac the distance of the cluster from the 
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galactic centre. We assumed a spherically symmetric density 
distribution for the Milky Way with a constant circular ve- 
locity of Vg ~ 220 km/sec. Cluster masses were calculated 
from the absolute luminosities of the clusters and an as- 
sumed V-band mass-to-light ratio of M/Ly = 2.0. The clus- 
ter data (projected half-light radii, total luminosities, galac- 
tocentric distances) was taken from the 200 3 version of the 
Globular Cluster database of iHarri^ (|l996), supplemented 
by ad ditional data for a few clusters from Bonatto Sz Bical 
l|2008h . We exclude Omega Cen, M54 and NGC 2419 from 
our analysis since these clusters might be stripped nuclei of 
dwarf galaxies rather than genuine globular clusters. 

Fig. [2] depicts the ratio of 3D half-mass radius to 
Jacobi radius r j as a function of galactocentric distance. It 
can be seen that the distribution of Galactic globular clus- 
ters is not uniform in this plane. First, clusters with large 
values of r^/rj > 0.5 are basically absent. Since such clus- 
ters would be subject to strong tidal forces and have con- 
sequently small dissolution times, if they ever existed they 
should be quickly destroyed and not be present any more af- 
ter 10 Gyr of evolution. The absence of such clusters in our 
distribution therefore serves as a sanity check of our method. 

Second, most clusters inside ~ 8 kpc exhibit a relatively 
broad distribution of r^/rj values between 0.02 < r^/rj < 
0.2 without any noticeable separation. Outside about 8 kpc, 
the Galactic globular clusters can be split into two groups, 
one group of clusters with r^/rj < 0.05 and a second group 
of clusters with 0.08 < r^/rj < 0.3. A KS test gives only an 
14% chance that clusters beyond 8 kpc follow a log-normal 
distribution in log r^/rj. In addition, the average mass of 
extended clusters in the outer Milky Way is significantly 
lower than the mass of the more compact clusters (see Fig. 
4), the mean mass of compact clusters is log Mc = 5.44±0.11 
while the extended group has a mean mass of only log Mc = 
4.38 ± 0.06. Both results indicate that two distinct groups 
of globular clusters exist in the Milky Way. The dividing 
line between both groups seems to be around r^/rj = 0.07 
and Table 1 lists the basic parameters of clusters having 
rfi/rj ratios smaller or larger than this value. Different orbits 
seem unlikely to be an explanation for this dichotomy since 
clusters of both groups are located within the same interval 
of galactocentric distances and, at least for those clusters 
with orbital information, the average ratios of Rac/Rperi 
are similar. 



The dashed line in Fig. [2] depicts the position which a 
Mc = 10* Mq with rh = 3 pc would have in the rh/r,j vs. 
Rgc plot at different galactocentric distances. It can be seen 
that clusters in the lower group outside 8 kpc and most clus- 
ters inside this radius fall onto this line. Most clusters in the 
lower group are therefore relatively massive and compact, 
something which can also be seen in Table 1. Clusters in the 
upper group strongly feel the tidal field of the Galaxy. It can 
be seen that this group is more diverse since it is made up 
of massive and extended as well as low-mass compact clus- 
ters. In the following, we will discuss possible reasons for the 
origin of both groups. 




Rgc [kpc] 

Figure 2. Ratio of half-mass radius to Jacobi radius rj as 
a function of galactocentric distance Rqc for Galactic globular 
clusters. It can be seen that clusters outside Rgc ~ 8 kpc fall 
into two distinct groups, clusters with rf^/rj < 0.05 and clusters 
with 0.07 < ri^/rj < 0.3. The dashed line depicts the position 
that a Mc = 10^ Mq with = 3 pc would have at different 
galactocentric distances. Compact clusters outside 8 kpc and most 
clusters inside this radius fall onto this line. 



3 DISCUSSION 

3.1 Correlation with old and younger halo 
membership 

The Galactic Globular cluster system consists of different 
subsystems. While bulge/disc GCs differ from halo clusters 
with respect to their metallicity ([Fe/H] ^ —0.8 and [Fe/H] 
< —0.8, respectively), the halo subsystem itself is made of 
clusters with more than one origin. It is traditionally split- 
ted up into two groups, referred to as the Old Halo and 
the Younger Halo (Zinn 1993, van den Bergh 1993, Mackey 
& Gilmore 2004), based on differences in horizontal branch 
(HB) morphology, age, kinematics, spatial distribution (see 
Parmentier et al. 2000, their Section 2, for a review). Because 
of their predominant location beyond the Solar Circle, YH 
clusters are assumed to have been accreted - along with the 
dwarf galaxies which used to host them - after the main 
body of the Galaxy was built up. Depending on how late 
they were accreted into the Galactic halo, their evolution- 
ary history may be different from what m-situ OH clusters 
have experienced in the Milky Way tidal field. The current 
accretion of the Sagittarius dwarf galaxy and of its small 
globular cluster system is the smoking gun of this process 
(Ibata, Gilmore & Irwin 1994). 

In Fig. O data points are symbol-coded to highlight 
these different cluster origins. Filled circles depict disc GCs 
([Fe/H] ^ —0.8), open triangles show GCs associated to the 
merging dwarf galaxy Sagittarius (Ter7, Arp2, Ter8, Pall2, 
NGC4147; see Da Costa & Armandroff 1995, Martinez- 
Delgado et al. 2002 and Bellazzini et al. 2003 for cluster 
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Table 1. Basic data for Galactic globular clusters with Rqc > 8 kpc belonging to both groups 



Name 




rh/rj 


ogMc 


logTRjj 


Rgc 




Name 






log Mc 


logTflH 


Rgc 


Rperi 




[pc] 








[yr] 


[kpc] 


[kpc] 




[pc] 




[Mq] 


[yr] 


[kpc] 


[kpc] 






Compact 










Tidally filling group 












2.67 





023 


5.60 


9.09 


9.4 


0.8 


NGC 288 


7.58 


0.093 


4.93 


9.50 


12.0 


5 3 


NdC 1 2fil 


4.77 





032 


5.36 


9.37 


18.2 




Pal 1 


2.87 


0.104 


3.22 


8.23 


17.0 




Pal 2 


7.17 





029 


5.44 


9.67 


35.4 




AM 1 


23.64 


0.114 


4.12 


9.92 


123.2 






2.44 





015 


5.57 


9.02 


16.7 


5.7 


Eridanus 


13.99 


0.070 


4.29 


9.65 


95.2 




NGC iqr)4 


4.00 





026 


5.38 


9.27 


18.8 


4.2 


Pyxis 


20.79 


0.150 


4.53 


10.00 


41.7 




i'l v,.i v_/ Zi^yo 


3.24 





038 


4.75 


8.87 


15.7 


1.9 


Pal 3 


23.73 


0.100 


4.51 


10.08 


95 9 


82.5 


ly V .1 V_/ ZOUO 


2.83 





016 


5.99 


9.29 


11.1 


2.6 


Pal 4 


22.87 


0.079 


4.64 


10.10 


111 8 




NOP '^901 


5.20 





062 


5.22 


9.37 


8.9 


9.0 


Rup 106 


9.04 


0.094 


4.77 


9.55 


18.5 




Nnr; 414.7 


3.22 





032 


4.70 


8.85 


21.3 


4.1 


NGC 5053 


22.26 


0.218 


4.92 


10.20 


16.9 




NOP 4'iqn 


6.13 





070 


5 17 


9.46 


10.1 


8.6 


AM 4 


4.87 


0.175 


2.87 


8.46 


25.5 




NOP ^(124 


7.66 





039 


5.71 


9.83 


18.3 


15.5 


NGC 5466 


13.87 


0.130 


5.02 


9.93 


16.2 


5.4 


Nnr; ^272 


4.52 





028 


5.81 


9.52 


12.2 


5.5 


IC 4499 


11.00 


0.094 


5.17 


9.84 


15.7 




NnC '^9Sfi 

IN VJV_> O^OU 


2.94 





026 


5.68 


9.19 


8.4 




Pal 5 


26.63 


0.395 


4.30 


10.07 


18.6 


6.1 


NGC 5634 


5.28 





033 


5.31 


9.42 


21.2 




Pal 14 


32.96 


0.233 


4.13 


10.14 


69.0 




NGC 5694 


4.44 





022 


5.36 


9.32 


29.1 




NGC 6101 


10.15 


0.124 


5.00 


9.72 


11.1 




NGC 5824 


4.47 





017 


5.77 


9.50 


25.8 




Pal 15 


20.93 


0.175 


4.43 


9.96 


37.9 




NGC 6205 


4.45 





037 


5.71 


9.47 


8.7 


5.0 


NGC 6426 


7.71 


0.084 


4.91 


9.50 


14.6 




NGC 6229 


4.36 





020 


5.45 


9.35 


29.7 




Ter 7 


8.73 


0.148 


4.25 


9.33 


16.0 




NGC 6341 


3.47 





031 


5.51 


9.23 


9.6 


1.4 


Arp 2 


21.19 


0.276 


4.35 


9.94 


21.4 




NGC 6779 


4.54 





053 


5.19 


9.27 


9.7 


0.9 


Ter 8 


10.08 


0.152 


4.25 


9.42 


19.1 




NGC 6864 


3.77 





023 


5.65 


9.34 


14.6 




Pal 12 


9.48 


0.193 


4.03 


9.29 


15.9 




NGC 6934 


3.65 





034 


5.22 


9.14 


12.8 


6.0 


Pal 13 


4.60 


0.083 


3.73 


8.71 


26.7 




NGC 6981 


5.80 





062 


5.05 


9.37 


12.9 




NGC 7492 


12.21 


0.124 


4.54 


9.66 


24.9 




NGC 7006 


6.12 





026 


5.31 


9.51 


38.8 


18.2 


IC 1257 


13.57 


0.149 


4.69 


9.79 


18.5 




NGC 7078 


4.23 





027 


5.90 


9.52 


10.4 


5.4 


BH 176 


4.84 


0.138 


3.97 


8.84 


10.2 




NGC 7089 


4.15 





028 


5.84 


9.48 


10.4 


6.4 


ESQ 280 


8.42 


0.156 


4.19 


9.28 


15.0 





membership). Plus-signs stand for GCs with no HB mor- 
phology index. A list of OH clusters (filled squares) is pro- 
vided in Parmentier & Grebel (2005, their Table 1). Other 
clusters are sorted in the YH group (open squares). Mackey 
& Gilmore (2004) emphasize that accreted clusters could 
also contribute a small fraction of the OH component. Based 
on either large core radius reminiscent of those observed for 
GCs in satellite galaxies (their fig. 16) or spatial motions 
more typical of YH objects (see also Dinescu et al. 1999), 
they identify 11 OH GCs which might have been accreted 
(NGC 6809, 6101, 7492, 5897 and Pal 15 in the first cat- 
egory and NGC 1904, 2298, 5024, 5904, 6205, 7089 in the 
second category). These ill-defined status clusters are shown 
as filled squares with open circles in Fig. [3] Clusters from 
Bonatto & Bica (2008) finally are shown as filled triangles. 
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It can be seen that there is no correlation among the 
classification of a cluster to either younger or old halo group 
and its rh/rj value. The compact cluster group contains 11 
younger halo clusters and 10 old halo clusters unsuspected 
of having been accreted. They represent each ~ 40 % of the 
total number of clusters (26) in the compact cluster group. 
The tidally filling cluster group contains 10 younger halo 
clusters and 5 old halo clusters unsuspected of having been 
accreted. The higher fraction of younger halo clusters, how- 
ever, is mostly driven by the 4 clusters with Rgc ~ 100 kpc. 
Considering the same radial extent as for the compact group, 
that is, Rgc = 8 — 50 kpc, younger and old halo clusters 
contribute similarly to the tidally filling group, with 6 and 
5 clusters, respectively, out of 22 clusters. Corresponding 



Figure 3. Ratio of r^/rj vs. galactocentric distance for Galac- 
tic globular clusters. Different symbols indicate to which sub- 
system a cluster belongs, old halo clusters are marked by filled 
squares, younger halo clusters by open squares, disc globular clus- 
ters by filled circles and clusters accreted from the Sagittarius 
dwarf galaxy by open triangles. There is no correlation between 
the classification of a globular cluster and its r^/rj value. 



number fractions (~ 25 %) agree with previous ones within 
the statistical uncertainties. Moreover, both compact and 
tidally filling groups are characterized by the same num- 
ber fraction of clusters accreted or suspected of having been 
accreted (i.e. younger halo clusters or Sagittarius clusters 
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or old halo clusters suspected of having been accreted - see 
above), namely, ~60%. We therefore conclude that the rh/rj 
dichotomy is not due to a different origin of the two cluster 
populations 



3.2 Origin of the compact cluster group 

Tab. [T] lists the basic parameters of clusters belonging to ei- 
ther group, including the 3D half-mass radius, current mass 
and galactocentric distance. It also shows the cur rent relax- 
ation time, calculated according to lSpitzerl (Il987l l: 



Trh = 0.138 



3/2 



VG <m> InO.llMc/ <m> 



(2) 



where <m>— O.4M0 is the average mass of stars and G 
the gravitational constant. It can be seen that the compact 
clusters mostly have very large relaxation times. The average 
relaxation time for a cluster in this group is about 2.8 Gyr, 
and nearly all clust ers have relax a tion t imes larger than 1 
Gyr. According to iGiirkan et al.l (120041 ) . it takes about 7 
to 10 initial half-mass relaxation times until star clusters 
with a narrow mass spectrum where the massive stars are 
about twice as massive as the average cluster star, which 
is typical for globular clusters, have gone into core collapse. 
The compact clusters should therefore still be mostly in their 
pre-core collapse phase and should not have started post- 
core collapse expansion. 

Tab. [T] also lists the perigalactic distances of 
the globul ar clusters as determined by Dinescu et al.l 

(|2006l ) 



(1 19991) 



Allen. Moreno fc Pichardol 
20071 ). 



and 

ICasetti-Dinescu et al.l (|2007| ). Although there are clus- 
ters which have perigalactic distances of less than 2 kpc, 
for the majority of the compact group clusters, the peri- 
galactic distances are within a factor of 3 of the current 
galactocentric distance. The estimates of r j would therefore 
decrease by no more than a factor of 2 if we used the 
perigalactic distance to calculate rj. Most compact clusters 
therefore have rn/rj < 0.1 also at perigalacticon and are 
at most moderately influenced by the Galactic tidal field. 
Hence, their small half-mass radii are likely not due to tidal 
stripping at perigalacticon but must have been the result of 
the formation process. 

We conclude that clusters in the compact group also 
formed very compact. A^'-body simulations show that the 
expansion factor due to gas expulsion is typically a factor 
2 to 3 for moderate star format ion efficiencies of 30% to 
40% l|Baumgardt fc Kroupa|[2007l ) and stellar evolution will 
increase this value by another factor 2 if mass is lost adia- 
batically. Hence the initial half-mass radius of clusters in the 
compact group must have been around 1 pc or less. Since 
most star clusters inside 8 kpc have half-mass radii very 
similar to compact group star clusters, it seems likely that 
most globular clusters in the Milky Way formed compact 
and with half-mass radii of 1 pc or less, which is compara- 
ble to the h alf-mass radius of em bedded star clusters in the 
Milky Way (|Lada fc Ladall2003l ). 



3.3 Origin of the tidally filling cluster group 

Fig. |4] depicts the position of inner clusters (green triangles) 
and of clusters with ru/rj < 0.07 (red crosses) and clusters 



with rh/rj > 0.07 (blue dots) in a half-mass radius vs. mass 
diagram. It can be seen that clusters with rh/rj < 0.07 
are mostly massive clusters with half-mass radii of a few 
pc, while clusters with rh/rj > 0.07 have larger radii and 
also smaller masses. Due to the smaller masses, clusters 
in the tidally filling group should on average be closer to 
dissolution. This is confirmed by observational data for a 
few clusters like Pal 5, which has very pronounced tidal 
tails and might be on i ts final orbit befor e dissolution 
dOdenkirchen et aLlliooH : iDehnen et al.ll2004l ). Clusters in 
the inner Milky Way also have smaller masses than compact 
outer clusters which might be due to stronger cluster disso- 
lution in t he inner Milky Way as a result of the stronger 
tidal field l|Vesperini fc Heggielll997l : iBaumgardt fc Makinol 
l2003l ). 

One way to explain the large radii of the tidally filling 
clusters would be that they also formed extended. Indeed, 
lElmegreenI (|2008l ) has recently discussed different modes of 
star formation and attributed the difference between star 
formation in bound clusters and loose groupings to a dif- 
ference in cloud pressure and different background tidal 
forces. This could explain why clusters with low densities 
are only found far away from the centers of major galax- 
ies or in dwarf galaxies. The fact that Milky Way globu- 
lar clusters are clearly separated in rn/rj is however more 
difficult to understand if cluster radii are set at formation 
time. An alternative viewpoint would be that the tidally 
filling clusters expanded from smaller radii, possible e.g. 
through post-collapse expansion dri ven by a p o pulati on of 
stellar binaries in the cluster core. iGoodmanI (|l984l ) and 
IBaumgardt. Hut fc Heggiel (|2002l ) (their Eq. 4) estimated 
that during post-core collapse expansion, the half-mass ra- 
dius of an isolated cluster satisfies 



rh{t) = r/io {t/tc 



\(2-l-f)/3 



(3) 



where Vho is the initial half-mass radius, tec the time of 
core collapse and ~ 0.1 a constant rela ted to the clus- 
ter mass loss. iGieles fc B aumgardtl l|2008l ) found that the 
above relation also holds for clusters in a tidal field as long 
as Th/rj < 0.05. For clusters with a narrow mass spectrum, 
core collapse happens after 7 to 10 initial half-mass relax- 
ation times (iGiirkan et al.ll2004l ). in which case the above 
relation would predict that expanding clusters should have 
relaxation times which are roughly 1/lOth of their current 
age, i.e. of order Trh ~ 10^ yrs. The majority of clusters 
in the tidally filling group however have relaxation times 
Trh > 3 • 10^ yrs, which is too large to be explained by 
binary driven expansion from small radii. Also the fact that 
tidally filling clusters have on average larger relaxation times 
than compact clusters argues against post-collapse expan- 
sion from smaller radii. 

Merritt ct al. (2004|) and iMackev etal] l|2007l ) have 
shown that stellar mass black holes, if present in sufficient 
numbers, can cause strong cluster expansio n. For clusters 
retaining all the black holes formed in them, iMackev et al.l 
l|2008t ) found that the core radius can reach values up to 8 pc 
after 10 Gyr of evolution and is almost as large as the half- 
mass radius. This value is large enough to explain the half- 
mass radii of a significant fraction of clusters in the tidally 
filling group (see Fig.|4]). Interestingly, in such a case clusters 
of the tidally filling group would have been the most com- 
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pact clusters initially such as to be able to retain their BHs. 
However, some clusters in the tidally-filling group have half- 
mass radii too large to be explained by BH-driven expansion 
and there is no significant difi^erence in metallicity between 
compact and tidally filling clusters, as might be expected 
if BH kick velocities depend on metallicity, which both ar- 
gue against BH driven expansion. Central intermediate-mass 
black holes can also act as an efficient heat source, but 
j udgin g from the results of lBaumgardt. Makino fc Ebisuzakil 
l|2004 ). the half- mass radii of most clusters in the tidally fill- 
ing group are too large to be explained by intermediate-mass 
black hole driven expansion. Strong expansion is also pos- 
sible by stellar evolution if star clusters are initially mass 
segregated since the fractional loss of potential energy can 
in such a case be much larger than the mass fraction lost by 
stellar evolution l|Vesperini et al.ll2009l ). 

The question of whether clusters in the tidally filling 
group were born compact and later expanded or already 
formed with the large half-mass radii we see today therefore 
remains open. If they formed with large half-mass radii, their 
initial relaxation times were also quite large and the clusters 
should be dynamically less evolved. In this case they would 
not be mass segregated, so measuring stellar mass functions 
at different radii might be one way to test the formation 
scenario. In this context it is interesting to note that Jordi 
et al. (2009) recently found that the stellar mass function of 
Pal 14, which is one of the clusters with the longest relax- 
ation time in our sample, differs from a Kroupa IMF inside 
the clusters half-mass radius. Unfortunately no information 
on the stellar mass function in the outer cluster parts is avail- 
able at the moment to test whether this is due to dynamical 
cluster evolution. 

We note that the tidally filling clusters, due to the 
large half-mass radii and galactocentric distances of many 
of them, are partly responsible f or driving the correlation 
betwe en th and Rgc found by iMackev fc van den Bergbl 
l|2005l ). A linear least-square fit gives a relation between 
half-mass radius and galactocentric distance log = 0.50 ■ 
log Rgc + 0.27 for all Galactic GCs. This relation fiattens 
to log rh = 0.25 • log Rgc + 0.38 if we exclude the tidally fill- 
ing clusters, indicating that globular clusters in the compact 
group formed with similar parameters nearly everywhere in 
the Galaxy. 



4 CONCLUSIONS 

We have studied the distribution of the ratio of half-mass 
radius to Jacobi radius r j for Galactic globular clusters 
and have shown that clusters with distances larger than ~ 
8 kpc fall into two distinct groups: One group of compact 
clusters with r^/rj < 0.05 and a group of more extended 
clusters with 0.08 < rn/rj < 0.30. Compact clusters are 
mainly massive clusters with half-mass radii of a few pc. 
The half-mass radius and density of the compact clusters 
in the outer halo seems not to be adjusted to the Jacobi 
radius, so they were probably also born compact with half- 
mass radii < 1 pc, comparable to the half-mass radii 
of embedded clusters and young open clusters in the Milky 
Way. Tidal radii derived from fitting King profiles to the 
surface density profiles of these clusters can be significantly 
smaller than their Jacoby radii since the Vh/rj ratios of these 
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Figure 4. Half-mass radius r;, vs. cluster mass for globular 
clusters with Rgc < 8 kpc (green triangles) and for clusters 
Rgc > 8 kpc that are weakly influenced by the Galactic tidal field 
('f'h/fJ < 0.07, red crosses) and for strongly tidally influenced 
clusters (0.07 < ffj/rj < 0.3, blue dots). Clusters weakly influ- 
enced by the tidal fleld are all massive and compact while strongly 
tidally influenced clusters have signiflcantly smaller masses. Inner 
clusters also have smaller masses on average, which might be a 
result of their stronger dissolution. Dashed lines show where clus- 
ters with a given half-mass relaxation time are located in this 
plot. 



clusters are smaller than what can be reached with any King 
profile. 

Some of the tidally filling clusters might also have 
formed compact and could have expanded later due to dy- 
namical heating by binary stars, stellar-mass black holes or 
intermediate-mass black holes, although it is unclear if this 
holds for all clusters in this group since about half of the 
tidally filling clusters have relaxation times of the order of 
a H ubble time or lar g er. 

iDa Costa et al.l (|2009l ) have recently found a bi- 
modality of the globular cluster size distribution in dwarf 
galaxies. The average radii of clusters in both of their groups 
agree quite well with the radii of Milky Way globular clus- 
ters in our compact and tidally filling group, showing that 
globular clusters forme d under similar c onditions in different 
galaxies. Furthermore, iPfalzneil l|2009l ) has recently shown 
that open clusters in the Milky Way evolve along two se- 
quences in the age vs. radius plane, one group of clusters 
starting compact with half-mass radii < 1 pc and reach- 
ing sizes of a few pc after 20 Myr of evolution and a second 
group of clusters starting with half-mass radii larger than a 
few pc and reaching ~ 20 pc after 20 Myr. The latter value 
agrees quite well with the sizes of most clusters in the tidally 
filling group. Galactic globular clusters therefore show the 
same dichotomy seen for globular clusters in dwarf galaxies 
and for young star clusters in the Milky Way. Extended star 
clusters appear therefore as an ubiquitous feature of star 
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cluster systems hosted by a variety of galaxies. It would be 
interesting to see how the extended globular clusters of the 
Milky Way relate to other star clusters with large rn val- 
ues, like the Faint Fuzz y star clusters in lenticular galaxies 
IILarsen fc Brodiell2000l ). the diffuse star clusters found by 
IPeng et alJ ( 2006f ) in early- type gala xies of the Virgo clust er 
and those hosted by dwarf galaxies JPa Costa et al.ll2009l ') . 

iKiipper. Kroupa fc Baumgardtl (|200^) found that ini- 
tially compact star clusters in a tidal field expand after core 
collapse until they reach a mass-dependent rn /rj value and 
then evolve along a common sequence towards dissolution. 
They dubbed the latter phase the main sequence evolution 
of star clusters. During the main sequence phase, the rn/rj 
values increase slowly with decreasing cluster mass. Extrap- 
olatin g from the results of iKiipper. Kroupa fc Baumgardtl 
l|2008l ) to Mc = 10^ Mq, we expect that globular clus- 
ters should have approximately Vh/rj « 0.1 when on the 
main sequence, which fits observed rn/rj of clusters in the 
tidally filling group rather well. It is therefore likely that 
part of the clusters in the tidally filling group, especially 
those with small relaxation times, have reached the main- 
sequence stage of their evolution and are evolving towards 
dissolution. As for those with long relaxation times, however, 
whether their large half-mass radius is an imprint of their 
formation process or a result of cluster expansion remains 
an open question. 

We finally note that it is possible that the extended clus- 
ters were initially much more numerous, since due to their 
large sizes, they are effectively destroyed by the Galactic 
tidal field, especially in the inner part of the Milky Way. 
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